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Abstract 

The recent observations of /o(980) in charmless -B-decays motivate 
further studies of scalar particle and glueball production in these pro- 
cesses. Amplitudes for charmless 2-body B decays involving the mem- 
bers of the scalar nonet are presented based on the symmetries of the 
dominant penguin contribution. Different scenarios for the lightest 
scalar nonet are investigated in view of the presently available data. 
We describe the evidence from £?-decays for /q(1500) with a flavour 
octet like mixing and the hints towards the members of the qq nonet 
of lowest mass. There is further support for the hypothesis of a broad 
++ glueball acting as coherent background especially in B — > KKK. 
The estimated B decay rates into gluonic mesons represent a sizable 
fraction of the theoretically derived decay rate for b — > sg. 
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1 Introduction 



A longstanding problem in QCD is the prediction of gluonic boundstates 
( "glueballs" ) and the lack of their doubtless identification. The lightest glue- 
ball is expected in the scalar channel with J PC = ++ . In this channel there 
is a series of established resonances but also various broad objects whose 
existence is not generally accepted. There is also no consensus about the 
lightest scalar qq nonet, neither about its members nor about the mixing 
between strange and nonstrange components. A central issue is the nature 
of / (980) which has been considered not only as standard qq meson but also 
as i^K-molecule or as qqqq state. 

There are now new experimental results on charmless B decays into scalar 
particles which provide additional information of high statistics. In this paper 
we discuss some recent results and their implications on scalar spectroscopy: 

1. The observation by the BELLE jl!2l.']l4j and BaBar collaborations jr>l(ij 

of charmless decays B — > Kh + h~ with h = tt,K with a significant peak 
related to / (980) and some less pronounced signals from other scalars. 

2. In the channel B — > KKK BELLE has also observed a broad enhance- 

ment (coherent "background" ) in the KK mass spectrum in the range 
1000-1700 MeV and perhaps beyond with spin J = and a smaller 
effect in inr around 1000 MeV. Similar results are shown by BaBar but 
there is no quantitative analysis yet. 

Remarkably, in B decays, as compared to decays of charmed D or J / 'ijj mesons 
the role of scalar particles is more pronounced because of the larger phase 
space and therefore reduced background from crossed channels and also be- 
cause of an apparent suppression of higher spin ( J = 2) states. In this paper 
we will discuss what can be concluded from the above observations and how 
further experimental studies can uncover the members of the light scalar 
nonet besides / (980) and possibly help to identify the gluonic interaction 
related to the ++ glueball. Some considerations concerning earlier data have 
been presented already in [7] . Other pertinent results concern the observation 
of /o(980) and gluonic mesons in gluon jets |7|8| . Not much theoretical work 
so far has been concerned with 5-decays into scalar particles. Some selected 
5-decays into flavoured scalar mesons have been derived from a factoriza- 
tion ansatz for the effective weak Hamiltonian including power corrections 
0(A/m b ) 0. 

The interest in charmless 5-decays with strangeness has been stimulated 
through the observation by CLEO |l()fllj of large inclusive and exclusive 
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decay rates B — > rj'X and B — ► f/lf, which have been confirmed by more 
recent measurements |12|13|l4j . These processes have been related to the 
decay b — > sg of the 5-quark which could be a source of mesons with large 
gluon affinity |15pi6|17|18j . In consequence, besides 77' also other gluonic 
states, in particular also scalar mesons or glueballs could be produced in a 
similar way. 

The total rate b — > sg has been calculated perturbatively in leading [19] 
and next-to-leading order [201 



„ /L , f(2-5)xl0- 3 inLO (for // = m 6 . . . m b /2) n . 
B(b^sg) = < . (1) 

I (5 ± 1) x 10 d m NLO 

The energetic massless gluon in this process could turn entirely into gluonic 
mesons by a nonperturbative transition after colour neutralization by a sec- 
ond gluon. Alternatively, colour neutralization through qq pairs is possible as 
well. This is to be distinguished from the short distance process b — > sqq with 
virtual intermediate gluon which has to be added to the CKM-suppressed 
decays b — > qi&Qz- These quark processes with s have been calculated and 
amount to branching fractions of ~ 2 x 1CT 3 each |21|22|2Uj . The question 
then arises which hadronic final states correspond to the decay b — > sg. 

We recall that a large gluonic penguin component has been suggested to 
play also an important if not dominant role in the explanation of the AI = ~ 
rule in K decay 

Next we outline the status of the phenomenological discussion concerning 
the light scalar spectroscopy. The Particle Data Group [21] lists below 1800 
MeV the following isoscalar particles 

/ Q (600)(or<7), / (980), / (1370), /„(1500), / (1710) (2) 

furthermore the isovectors a (980), a (1450) and the strange ^(1430), pos- 
sibly there are k(850) and _K"g(1950). The broad states <j,k are still contro- 
versial, also not much reliable information is available about /q(1370). The 
states listed above should be related to a scalar qq nonet and, possibly, a 
glueball. 

Quantitative results on glueballs are derived today from the QCD lattice 
calculations or QCD sum rules, both agree that the lightest glueball has 
quantum numbers J PC = ++ . 

Lattice calculations in quenched approximation |25|2fi|27|28] (without 
light sea quark- ant iquark pairs) suggest the lightest glueball to have a mass 
in the range 1400-1800 MeV [2*9*] . Results from unquenched calculations still 
suffer from systematic effects, the large quark masses of the order of the 



2 



strange quark mass and large lattice spacings. Typically, present results on 
the glueball mass are about 20% lower than the quenched results |30|31| . 

Results on glueballs have also been obtained from QCD sum rules. Re- 
cent calculations [SHj for the ++ glueball yield a mass consistent with the 
quenched lattice result but in addition require a gluonic state near 1000 MeV. 
Similar results with a low glueball mass around 1000 MeV are obtained also 
in other calculations @2]. On the other hand, it has been argued [41J that 
the sum rules can also be saturated by a single glueball state with mass 
1250 ± 200 MeV. 

In conclusion, there is agreement in the QCD based calculations on the 
existence of a ++ glueball but the mass and width of the lightest state is 
not yet certain and phenomenological searches should allow a mass range of 
about 1000-1800 MeV. 

In the interpretation of the phenomenological results one can identify two 
major different directions of thought; they differ in the role of /b(980) which 
belongs either to a nonet with particles of higher mass or with particles of 
lower mass. 

Route A: Scalar nonet with / (980) and heavier particles 
In a previous study [S21 we have performed a detailed phenomenological 
analysis of production and decay of low mass scalar mesons, which led us to 
identify the lightest qq scalar nonet with the states 

/o(980), a (980), K*(1430), / (1500) (3) 

with large flavour mixing, just as in the pseudoscalar nonet, i.e. with flavour 
components (uu, aid, ss) approximately given by 

rf, / (980) <-> (1,1,2)/V6, V , / (1500) <- (1,1,-1)/^ (4) 

close to the flavour singlet or octet respectively and with the parity partners 
rf and / (980). 

This scalar nonet separates to a good approximation into the singlet 
/o(980) and the octet a (980), #5(1430), / (1500). The octet, within this 
approximation, fulfills the Gell Mann-Okubo mass formula and is also con- 
sistent with a chiral model with general QCD interaction of E fields. In 
elastic tttt scattering and in other channels the cross section shows a broad 
"background" which extends from below 1000 MeV up to about 1600 MeV or 
higher, with dips from negative interference with narrow states /o(980) and 
/o(1500) ("red dragon"). The PDG has been listing here the two states left 
over so far / (600) and / (1370). We have interpreted this "background" as 
signals from a single broad object /o(1000) centered around 1000 MeV with 
a large width of 500-1000 MeV (according to the T matrix pole parameters 
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|33|34| ). This "left over state" was taken as ++ glueball as its production 
characteristics followed largely the expectations for a gluonic meson (except 
for J /if) — > 77T7r). Subsequent studies of various decay rates and relative 
amplitude phases involving / (980) and /o(1500) have confirmed this view 
|35|36| . sometimes better experimental data would be desirable. The light 
scalar objects a, k are not considered as individual physical particles to be 
classified into a nonet. The heavier particles a (1450), / (1710) could be 
members of a second nonet, possibly with if* (1950) and /o(2020). 

There are other approaches which agree on /o(980) as being the lightest 
particle in the nonet. A scheme similar to ours above for the qq nonet with 
/o(980) and / (1500) and mixing like (@J) has been proposed [37] based on 
a quark model with instanton interactions. It prefers though do (1450) over 
a (980) as nonet partner of / (980), but there is no prediction on a glueball. 
The broad enhancements in tttt and also An spectra have been interpreted as 
a consequence of p-exchange in two-body scattering amplitudes |3*H] . 

A broad glueball (width about 2 GeV) is found in the mass range 1200- 
1600 MeV from overlapping f states in a if matrix analysis of a variety of 
reactions [12] • The glueball nature of this state has been inferred from its 
"flavour blind couplings" . Two qq multiplets emerge in the mass range below 
1900 MeV including /o(1300) and again /o(980) as the lightest scalar which 
now behaves like a flavour octet. 

Using results from QCD sum rules a scheme for light scalars has been 
proposed PH] where a broad isoscalar "a" around 1000 MeV and the narrow 
/o(980) are both mixed in equal parts from the broad glueball and a light qq 
scalar. 

Despite various differences in detail, especially the role of the glueball, all 
these schemes have /o(980) as the lightest member of the nonet. 
Route B: Scalar nonet with / (980) and lighter particles 
In a second line of thought there is one qq nonet at heavier mass and, in 
general, a second nonet qq or qqqq at lower mass including /o(980). The 
heavier multiplet includes ^(1430), the isovector candidate in the same mass 
region is a (1450). In the isoscalar channel one observes / (1370), / (1500) 
and /o(1710). The scalar glueball is assumed with mass around 1600 MeV 
as suggested by quenched lattice calculations. Then the glueball and two 
isoscalar members of the nonet can mix and generate the three observed /o 
states above. Several such mixing schemes have been proposed giving either 
a larger gluonic component to / (1500) [13] or to / (1710) jUj. 

The lower mass scalars are now left over. An interesting possibility is the 
existence of a light nonet including 

/o(980), a (980), k(850), a(600) (5) 
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(alternatively also if* (1430)). This nonet could be of conventional qq type 
|45|46|47|48| or built from qqqq |49|50|5lj . Mesons of this light nonet can 
also be related to poles in meson meson scattering amplitudes constructed 
using chiral symmetry and unitarity Accepting the existence of a and 
k particles one could have two nonets below 1800 MeV (for reviews, see also 
|53l54l55j ). 

Finally, we should comment on the problem of the a, k poles whose inter- 
pretation as particles remains controversial |56|57j . The evidence has been 
studied in detail in tttt scattering within a class of parametrizations which 
respect chiral symmetry and unitarity The small scattering length and 
the rapid increase of the isoscalar S wave phase shifts with energy required by 
unitarity yields a pole in the amplitude related to a. Whether the pole in this 
parametrization represents a physical particle is not immediately obvious. 

In particular, it has been pointed out jSHj that final state irn interactions 
calculated within the non-linear sigma model up to two loops "mock up" 
such a particle: the "<r" particle is nothing but a convenient parametrization 
of the correlated two pion exchange. Here we wish to argue that - if these 
objects are really physical particles - the decay amplitudes should respect 
the associated flavour symmetry relations which we discuss below. 

The aim of our paper is the identification of the lightest scalar nonet 
and possibly the light scalar glueball. The first step in the next Section is a 
classification of amplitudes in charmless B decays. We describe an approxi- 
mation for decays into two pseudoscalar (PP) or into a pseudoscalar and a 
vector particle (PV) employing U(3) flavour symmetry and compare with the 
recent experimental results; gluonic transitions play an important role. In 
Section 3 we turn to the discussion of scalar particles (S) in final states (PS) 
and (VS) which is the main interest of the paper and we present the corre- 
sponding decay amplitudes. For illustration we compare our formulae with 
the first experimental data for two different scenarios of the scalar sector and 
present further predictions. One possible interpretation of the data involves 
a broad scalar glueball; in Section 4 we estimate the total gluonic contribu- 
tions to B decays and compare with perturbative calculations. Conclusions 
follow in Section 5. The detailed study of charmless B decays should help 
distinguishing the various possibilities outlined above. 

2 Charmless B decays with K and if* (890) 

Our aim is the description of charmless B decays into scalar particles (S) 
together with particles from pseudoscalar (P) and vector (V) multiplets, B — > 
PS, B — > VS, including especially the recently observed scalar / (980). In 
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Figure 1: Short distance processes for charmless b decays (a) QCD penguin 
diagram with qq pair production through intermediate gluon (U represents 
u, c, t quarks) and (b) with spectator annihilation; (c) penguin diagram with 
hard gluon; (d) electroweak (CKM suppressed) tree diagram for qq pair pro- 
duction and (e) spectator annihilation. In the present analysis we keep only 
(a) and (c). 

view of the yet incomplete understanding of the scalar sector and the scarce 
data we are interested in the leading approximation which describes the main 
effects. We begin by reconsidering the well studied decays B — > PP, B — > 
VP, especially B —>■ Krf, K*rj', together with the other final states related 
by U (3) symmetry. Subsequently we extend these considerations towards the 
inclusion of scalar particles. 



2.1 Approximation for two-body decays 

The large branching fraction B — > Kr\' confirms the special role of rj' in 
these decays and it has been related |15pifipi7fT£] to the gluon affinity of 
rj', especially through the QCD axial anomaly which affects only the flavour 
singlet component. In the factorization approach many 2-body decays can 
be reproduced well but not B — > Krj' |6U|61j . Also, it appears difficult to 
explain the Krj' rate entirely by quark final states and the QCD anomaly 
within a perturbative framework [62 , a factor 2 remains unexplained. An 
improvement is possible by inclusion of QCD radiative corrections 63J but 
with considerable uncertainties. 
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Alternatively, one may introduce a phenomenological flavour singlet am- 
plitude which allows also for non-perturbative effects [TR. This amplitude is 
added to the dominant penguin amplitudes, the small tree amplitudes and 
electroweak penguins. Different decays are related by flavour £7(3) symmetry. 
Recent applications |64 66J of this scheme to 2-body B decays with strange 
and nonstrange pseudoscalar and vector particles yield a good overall agree- 
ment with the data in terms of a few phenomenological input amplitudes. 

Here we discuss first the 2-body B decays with K and K* in this way 
|(i4ffi(ij to understand the pattern of the observed rates and then extend the 
analysis to the scalar sector. For this purpose we restrict ourselves to a simple 
approximation. In the description of the short distance interaction we keep 
only the dominant QCD penguin amplitudes T q for b — > qqs, q = u,d,s with 
intermediate virtual gluon and T u = T d = T s as well as the penguin amplitude 
with hard gluon radiation (see Fig. Hfe.,c) while we neglect the electroweak 
tree diagrams which are CKM suppressed in \V US \ = 0.22 (Fig. HJi) as well 
as the spectator annihilation processes of any kind (Fig. |TJ),e) because of the 
large 6-quark mass. 

Concerning the suppression of the tree diagrams we refer to the calcu- 
lations |22|20| of decay rates at the quark level including penguin and tree 
amplitudes. The decays of interest to us are found with relative fractions [22J 
22% (b — ► uus), 18% {b — > dds) and 15% (b — > sss) of all charmless hadronic 
B decays (a total of (55 ± 5%), the remaining ~ 45% correspond to b — > sg). 
So the nonleading weak decay amplitudes which we neglect here modify the 
leading result for these rates from penguin amplitudes by about ±20%; they 
become essential, if CP violating effects are investigated. 

Next we discuss the hadronic two-body decays B — > xy with particles x, y 
belonging to U(3) multiplets X and Y which are related to the penguin dia- 
grams Fig. ^i,c kept in the present approximation. The b — > sqq amplitudes 
in Fig.GJi lead to hadrons either by connection of the spectator quark (g s ) in 
the B meson with the produced quark q or with the s quark as shown in Figs. 
Efe.,b. The second process is colour suppressed as the qq pair produced in a 
colour octet state cannot recombine directly into a hadron; we will neglect 
this contribution in the following. The hadronic penguin amplitudes p q xy for 
the decay B — > x(sq)y(qq s ) into particles from multiplets X and Y in Fig. [2^. 
are then assumed proportional within the given multiplet to the amplitudes 
T q in Fig. UK 

Ply = AX sq A qq 3 PXY, PXY = k XY T , T = T u = T d = T s (6) 

where AjL denotes the flavour coupling x — > qq' and hxy a hadronization 
constant. Likewise there are amplitudes with x and y exchanged. 
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(a) (b) 




Y Y 

(c) (d) (e) 

Figure 2: Diagrams for charmless hadronic two-body decays B — > xy (like- 
wise B — > yx) into hadrons x, y from flavour multiplets X, K related to QCD 
penguin amplitudes of Figs, la-c: with b — > sqq followed by (a) colour 
favoured and (b) colour suppressed hadron formation; with b — > sg followed 
by (c) gluonic colour neutralization and (d) qq colour neutralization or (e) 
glueball production. The dashed lines indicate quarks and gluons with typi- 
cally soft interactions. 

In addition, there is the contribution from the b — > sg penguin in Fig. 
which leads to hadronic final states either by non-perturbative soft colour 
octet (g, gg, . . .) or by colour triplet (qq) neutralization of the hard gluon as 
in Fig. Efc,d. 2 The amplitude in Fig. Efc contributes to the production of a 
meson y with flavour singlet qq component (such as the rf) which we write 
as 

s xy = A x S g a (A y uQ + A y di + A* s )sxy, with s X y = IxyPxy- (7) 

We take into account this contribution only for the production of scalar and 
pseudoscalar particles and neglect it for vector particles (uj,(f>). This can be 
justified by the smaller perturbative transition amplitude with exchange of 
at least three gluons in the second compared to two gluons in the first 

one; this difference is also considered responsible for the large flavour mixing 
in the pseudoscalar and the small mixing in the vector meson nonet [HH] . 
The smallness of this flavour singlet contribution has also been found in the 

2 For a further discussion of these two mechanisms and observable consequences, see, 
for example, Ref. [5]. 



8 



(a) (b) (c) 

Figure 3: Two-body decay B~ — > fT - ?]' with three amplitudes as in Fig. 
Ek,c: (a) amplitude p\_- , with s — > K~ , (b) exchange amplitude with 
s — > Vj and (c) amplitude sk- v ' for gluonic production of 77'. 

phenomenological analysis [66]. Furthermore, glueball production is possible 
through this process (Fig. Ek). The diagram in Fig. Eli has the same quark 
structure as the one in Fig. and therefore it is not kept independently. 

If the two particles belong to two different multiplets then p q yx with both 
decay particles interchanged, i.e. with s — > y, q s — > x, is written as P'xyPxy 
In case of -B-decays into an isoscalar meson with mixed strange and non- 
strange quark components (77,77', . . .) both amplitudes contribute and inter- 
fere which yields the full amplitude p q xy + /3xyPI v where the second term 
with Pxy = (—^-) L P'xy refers to the 2-particle state with reflected momenta 
(p — > —p) for orbital angular momentum L |67j . In particular, for B — > VP 
decays there is a relative (— ) sign for the interchange amplitude. 

The B branching ratios are then computed from the superposition of 
up to three amplitudes depending on the quark structure of x and y: the 
amplitude for decay B — > x(sq)y(qq s ), the exchange amplitude for decay 
B — ► y(sq')x(q'q s ) and the gluonic amplitude for decay into an isoscalar 
meson B — > x(sq s ) ?/(isoscalar qq). 

B(B -> xy) 

= \Pxy Pyx "I" ^ 

H Wl^A T * + 4^*/^ + A", Yl A l^vT\ 2 . 

q" 

with produced quarks q, q', q" . The decay B~ — > ^"77' where all three ampli- 
tudes contribute is shown in Fig. 01 For the nonets of pseudoscalar, vector 
and scalar mesons considered here we list the relevant parameters and some 
notations in Tabled 
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Table 1: Parameters for decays B — > xy (s — > x, q s — > y) for pseudoscalar 
(P), vector (V) and scalar (S) particles and some special choices in the present 



analysis 


XY 


penguin 






flavour singlet 


special choices 


PP 


Ppp 






spp = IppPpp 




VP 


Pvp 


Ppv = 


PvpPvp 


s V p = 'JvpPvp 


Pvp = —1, Ivp = Ipp 


PS 


Pps 


Psp = 


PpsPps 


sps = IpsPps 


Pps = ±1 










ssp = IspPps 




vs 


Pvs 


Psv = 


PvsPvs 


svs = IvsPvs 


Pvs = ~Pps, Ivs = Ips 



2.2 Decays into pseudoscalar and vector particles 

We apply this simple approximation based on the dominance of penguin 
diagrams first to the well studied B — > PP and B — > VP decays. Non- 
leading diagrams are considerd in [54f6"6] which are important for the study of 
CP violation effects (absent in our approximation). The relative magnitudes 
of the tree amplitudes are found there of the order of ~20%. 

Without loss of generality we can define pxy as real in our approximation, 
whereas Pxy and 7xy could be complex in general. However, we assume in 
the present application also Pxy and ^xy to be real, to begin with. Further- 
more, we assume equal recombination probability for the exchanged particles, 
i.e. \Pxy\ = 1 or Pxy = ±(— 1) L - These simplifications could be relaxed if re- 
quired by the data. The quark mixing parameters in the pseudoscalar sector 
are taken from (j3J), i.e. rj = (uu + dd — ss)/v^3 and rf = (uu + dd + 2ss)/\^6. 

Results for B — ► PP and B — ► VP decays are given in Table 121 the 
amplitudes entering Eq. (jHJ) in col. 2 and explicitly in col. 3 for T q — 1 
in units of ppp and pvp with parameters 7pp and jvp and Pvp = — 1. 
Repeating the calculation with Pvp = +1 would exchange the roles of K* + r] 
and K* + i]' and can therefore be excluded by comparing with the data. The 
B + decay rates are obtained by adjusting in each sector the normalization 
of |l?xy| 2 and by multiplying the amplitude squared in col. 3 with |pxy| 2 - 
The predictions for B° follow by multiplying |pxy| 2 with the ratio Tpo/Tp+ = 
0.921 |24j . We compare with branching ratios updated recently for B — > PP 
|S3! and B -> VP [SS] in col. 5. 

In a first approximation we adjust only the two penguin amplitudes ppp 
and pvp = appp from two rates (K°ti + , K*°tt + ) and set 7xy = 0. Then 
the overall pattern of the data is reproduced, except for the rates Krj' and 
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Table 2: Branching ratios for B + and B° decays into pseudoscalar (P) and 
vector (V) particles (cols. 4,5) and amplitudes in Eq. (JEJ) (col. 2), ^ypp, 7 yp 
and Pvp for gluonic and interchange processes; col. 3: T q ,a = Pvp/ppp set 
to 1; col. 4: a = 0.661, 7pp = j VP = 0.439, VP = -1, \p PP \ 2 = 20.6 x 10" 6 . 



B — > PP 


amplitudes 


T q = 1 


B th [10- 6 } 


#exp[10- 6 ] 


K°7l + 


T d 


1 


input ppp 


20.6 ±1.3 


K+n° 


J-T 

./o u 
V z 


i 

V ^ 


10.3 


12.8 ± 1.1 


K+TJ 


Mt u -t s + 1pp t) 

V o 


73 7pp 


1.3 


3.1 ±0.7 


K + 7]' 


±{T u + 2T s + A lPP T) 


^(3 + 4 7P p) 


input 7 pp 


77.6 ±4.6 




T u 


1 


19.0 


18.2 ±0.8 


K°n° 


V ^ 


i 

V ^ 


9.5 


11.2 ± 1.4 


K°T] 


^ r (T d -T s + lPP T) 


7s 7pp 


1.2 


< 4.6 


K°rf 


MT d + 2T s + A lPP T) 

v o 


75(3 + 4 7 pp) 


71.5 


60.6 ±7.0 


B ^VP 




a = 1, /?vp = —1 






K*°n+ 


ocTd 


1 


input a 


9.0 ± 1.4 


K* + 7l° 


a rp 
y/2 U 


i 


4.5 


< 31 


K* + T] 


-%(T U -p VP T s +7ypT) 


4=(2 + 7 yp) 


17.8 


25.9 ±3.4 


K* + T]' 


-%(T u + 2(3 vp T s + 4 1vp T) 


-4(-l +4 7vP ) 


0.9 


< 12 


p+K° 


aflvpTd 


-1 


9.0 


< 48 




-t= Pv pTI, 


1_ 

\/2 


4.5 


4.1 ±0.8 




V2' 


1 


4.5 


5.4 ±0.8 


<j)K + 




1 


9.0 


9.0 ±0.9 


K* + n- 




1 


8.3 


15.3 ±3.8 


K*%° 




1 


4.2 


0.4 ± 1.8 


K*°T] 


^(T d -A, P T s ± 7yP T) 


^(2 + 7vp) 


16.4 


17.8 ±2.0 


K*°7]' 


^(T d + 2^pT s + 4 7y pT) 


^(-1 + 4 7FP ) 


0.8 


< 6.4 


p-K + 


a/3ypT u 


-1 


8.3 


9.0 ±2.3 


p°K° 




i 


4.2 


< 12.4 


tuK° 


^PvpT d 


i 


4.2 


5.2 ± 1.1 


(j)K° 


aT s 


1 


8.3 


7.8 ±1.1 
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K*7] which are significantly too large by a factor of up to ~ 3 {K + r]'\ 30.9, 
K°rj': 28.5, K* + r]: 12.0, K*°r]\ 11.1). This conclusion follows here from 
flavour symmetry. The neglect of nonleading short distance terms, estimated 
to amount about 20% cannot account for the discrepency, which is then 
attributed to the additional flavour singlet amplitudes ^ PP ppp and JvpPvp- 

Predictions for real 7p = 0.439, determined from K + i]' are shown in col. 
4, in the VP sector we chose for simplicity 7yp = 7pp (see also the similar 
results in [HI])- At this level of approximation, with expected accuracy of 
~ 20%, there are no major discrepencies encountered. 3 

Whereas we do not intend to go beyond the present approximations, we 
may estimate the effect of chosing 7pp complex. Then we find from Kt] and 
Krf rates (Table 7pp ~ 0.67 exp(iip PP ), ip PP ~ ±67°. Using either real 
or complex 7pp we obtain an estimate of the gluonic part of the Krj', Krj 
production rate 

B(B + ^K + V ', K + V )\ gluonic = 3 \ lPP p PP \ 2 ~ (12... 28) x lO" 6 . (9) 

We conclude that the main effects are reproduced in each sector by two 
parameters, the penguin amplitude pxy and the flavour singlet parameter 
7xy, furthermore we have chosen f3y P = —1. 

Of particular importance for our further discussion are the large and small 
rates for Krj' and Krj respectively and the abundancies of K*i]' and K*rj in 
reversed order which the model explains after choosing f3yp = —1, so far 
consistent with the data. This is a consequence of the different sign of the 
exchanged amplitudes in the PP and VP multiplets |b'7|64| . a feature also 
present in other analyses |62|63j . 



3 The scalar nonet in B decays 

3.1 Branching ratios: expectations and observations 

We now turn to the main part of our paper, the study of scalar particle 
production in B-decays. This is motivated by the remarkably strong signal 
observed for the scalar meson /o(980) by the BELLE pQ (recent update [Sj) 
and BaBar Collaborations jS] in the decay B + — > K + -k + -k^ where almost one 
half of the total rate above background falls into this sub-channel 

B(B + ^K + f ( 980); f -> tttt) w 15 x 10" 6 , (10) 

3 There are two predictions with a ~2.5<r deviation from data (K + rj and K* + if)\ a better 
result could be obtained by fitting the three parameters to all branching ratios instead of 
three only. Alternatively, one may relax the condition for ^pp to be real or the condition 
Ivp = Ipp- 
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Table 3: Dominant contributions for B decays into scalar (S) + pseudoscalar (P) or 
vector (V) particles: penguin amplitudes pxy (normalized to 1 in each sector as in 
col. 3 of Table El), exchange and gluonic amplitudes /3ps,/3vs and •yps,7sp,lvs resp. 
with scalar mixing angle tps] in brackets results for sin cps = (ps ~ <Pp)', cols. 

3,6: upper sign for B°, lower sign B + . 



B ^ 


B+ -> 


normalization to 


B ^ 


B+ -> 


normalization to 


P + S 


P + S 


Pps 


V + S 


V + S 


Pvs 


K + a~ 


K°a+ 


1 


K* + a- 


K*°a+ 


1 


K°a° 


K + a° 




K*°a° 


K*+a° 




K°fo 


K + h 


^(1 + 2'jps) sirups 
+(Pps + Ips) cos ip s 
[^(l + 2/5 P5 + 4 7 p 5 )] 


K*°f 


K* + fo 


^(1 + 27y 5 ) siny2 5 

+{Pvs + Ivs) COS ip S 
[j E (l + 2p vs + A lvs )] 




K + fo 


^(1 + 2-fp S ) cos(f S 

-(Pps + Ips) simps 
{-L(l - (3 PS + lPS )] 






75(1 + 27^/5) cos^s 
-{Pvs + Jvs) sin Lps 




n+K? c 


Pps 




P + K? c 


Pvs 


n o K: o 






P°Kf c 


p'Kt 


T^Pvs 






^(-1 + Pps + 75p) 


"Kf c 










^(2 + Pps + ^sp) 






1 



for more details, see Table |U below. This large fraction for f^ + /o(980) is 
comparable to pseudoscalar decays and 3 times larger than K + p°. Such a 
large rate could be taken as a first hint at the gluonic affinity of this meson as 
well, but it is clear that a more definitive answer requires an analysis similar 
to the one with 7/ for scalar (S) particles as well on the basis of their flavour 
classification. Here we assume /o(980) to be a member of a U(3) nonet. 

In Table El we have written down the amplitudes for the decays B — > PS 
and B — > VS in the approximation as above, keeping only the QCD penguin 
amplitudes (pps,Psp,Pvs,Psv) and gluonic amplitudes (sps, ssp, sys) for 
the respective multiplets, see Table 1. We denote the states of the scalar 
nonet by a , fo, K* c , and /q, where the isoscalar states are mixed from 
strange and non-strange components as 

f = nn sin ip s + ss cos ips, f q = nn cos ip s — ss simp s (11) 
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with nn = (uu ± dd)/y/2 and mixing angle cp$- For a given nonet of scalar 
states Table 01 predicts the corresponding pattern of decay rates in terms of 
these parameters. 

For the PS decays the parameters are: the normalisation of the penguin 
amplitude pps and ^ps, Jsp, Pps and for VS decays the normalization 
Pvs an d 1vs,Pvs- According to our experience with the PP sector we may 
assume in the beginning real parameters (3, 7 and an equal recombination 
probability for exchanged processes and therefore restrict to l/Sps] 2 = 1. 
Furthermore, we assume, as in case of pseudoscalars, that the replacement of 
one pseudoscalar by one vector meson keeps the parameters 7, (3' unaltered 
and therefore we choose 

Pps = ±1, Ivs = Ips, Pvs = -Pps, ( 12 ) 

where the opposite sign for Pvs comes from the spin factor (— 1) L . 

Next we ask whether besides the very clear signal of /q(980) there is 
any evidence for production of other scalars in the data from BELLE (Ref. 
PQ and the recent update, still preliminary, with higher statistics [3 ) and 
BABAR 0: 

Table 4: B decays into scalars measured by BELLE [2j and BABAR [Sj 



showing statistical and systematic as well as model errors). 

B(B+ Z K+/ (980)) (10.3 ±l.l™J)xlO- 6 [3j 

xB{f -> 7T+7T-) (9.2 ±1.2 ±0.6ll;| ±1.6) x !0~ 6 [5J 

B{B + -> A- * (1430)tt + ) (25.0 ± 1.6±|f±?;g) x lO" 6 Fit C /I 

xB(K*° -> K+tt-) (6.00 ± 0.84±8;|i±8:i|) x 10^ 6 Fit C /II 

(25.1 ± 2.0 ± 2.9±g:| d= 4.9) x lO" 6 [5J 

B(B + -> ir + /o(1500)) (18.5 ± 0.5) x lO" 6 Fit B /l 

xB(f ^K + K-) (1.3 ±0.2) x 10- 6 Fit B /II 



K*(U30) 

A higher mass Kq has been seen by BELLE and by BaBar (BaBar quotes 
"higher K*" which includes higher spin states in this mass range) with rates 
reproduced in Table EJ Note the considerable difference in the two solutions 
by BELLE. 

/o(1500) and broad "background" in KK and tttt 

In the first publications by BELLE 1 2j there was a broad low mass en- 
hancement in the K + K~ and K°K° mass spectra without any further struc- 
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ture. The recent preliminary results t 3], show a qualitatively different pic- 
ture: there is a sharper peak at 1500 MeV in the K + K~ mass spectrum, 
whose mass and width agree with / (1500), above a broad bump ("back- 
ground"). In the 7r + 7r~ mass spectrum a broad "background" can be seen as 
well, however, there is no comparable peak at 1500 MeV which at first makes 
the identification with / (1500) difficult in view of the known branching ra- 
tio B(/ o (1500) -> K K)/B(f o (lh00) -> tttt) = 0.241 ± 0.028 [21] based on 
the measurements 68 69j . We will argue below that the identification with 
/o(1500) and its branching ratio is possible if the constructive and destruc- 
tive interferences respectively with the "background" are taken into account. 
In the K + K~ channel BELLE quotes for the resonance rate in two solutions 
Bo/ 1 and Bo/ II the quite different fractions 60.8% or 4.4% of the charm- 
less K + K + K~ final state from which we derive (with the errors from the 
fractions) the numbers in Table 0] 
/o(1370) 

Belle |3] also quotes the fractions of a small enhancement near 1300 MeV in 

7r + 7r~ as /x(1300) which could be a signal from / (1370). 

k(850) 

The existence of this state with parameters determined by the E791 collabo- 
ration [ZD! (Mass 797 MeV, width 410 MeV) could not be verified by BELLE 
[3]. If a resonance is fitted it would have the much larger width of 2.27 GeV. 
The best description of the decay B + — ► K + n + it~ by BELLE includes a 
broad background in Kir with constant phase. 
or(600) 

One may ask also whether there is any evidence for the <r(600) particle which 
is seen in some experiments as peak near the nir threshold. The recent 
BELLE data jHj do not show any peak below the p meson above background 
whereas the statistics in BaBar is too low to make definite statements. It 
would be interesting to have some limits or estimates of the B branching 
fractions into these hypothetical mesons a and k to be compared with the 
/o(980) branching fraction. 
a (980) 

An important state in scalar spectroscopy is the a (980) which directly mea- 
sures the penguin amplitude within the scheme of Table H3 A possible decay 
mode is a — > KK, which yields a peak just above KK threshold. It is 
interesting to note that the decay B° -> K + K~K° studied by BELLE |2j 
clearly shows such a threshold enhancement near 1000 MeV in K + K° and 
K~K°, which could be due to a (980) decay. However, there is a large back- 
ground in this region from non-5 decays and there is no claim for observation 
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of ao (980) by the BELLE collaboration. 4 A determination through the de- 
cay mode a (980) — > r/n would be interesting, also for the heavier do (1450). 
Such studies would decide whether the nonet partner of /o(980) is ao(980) or 
o (1450). 

The above list suggests the observation of scalar states heavier than 
/o(980) whereas there is no clear evidence yet for the observation of the 
states of lower mass. This is along route A which we discuss next. 

3.2 Comparison with heavy nonet scheme 

In this scheme /o(980) is the lightest member of the nonet. Specifically, in 
our version [32J the members are as in © and in addition there is a broad 
glueball. The scheme jHZ| is similar, with preference for ao(1450) but without 
glueball. In both schemes /o(980) and /o(1500) are mixed as given by (jlj). 

We begin with a discussion of the production of / (1500). The peak in the 
K + K~ mass spectrum and a lack of a signal in 7r + 7r~ finds a natural interpre- 
tation in the quark structure (J3J) with the negative sign between nonstrange 
and strange components |H2|H7| together with the glueball interpretation of 
the background [32] . This negative relative sign transfers into the decay am- 
plitudes of / Q — ► KK on one side and f — ► tttt or f — > r\r\ on the other 
side, whereas for the glueball the decay amplitudes have the same sign for 
all pseudoscalar pairs. Therefore, if the interference / (1500) - background 
is constructive for KK we expect it to be destructive in tttt. We recall that 
such a difference has been observed already in the channels tttt — > KK and 
tttt — > r\r\ in their interference with the broad background according to our 
earlier analysis [32*] . 

In order to check the validity of this proposal quantitatively, for the given 
decay ratio of /o(1500) into KK and tttt we compare the data with a simple 
parametrisation of the decay amplitudes in terms of /o(980), /o(1500) and the 
broad background which we interpret as glueball (gb). The new feature of this 
fit in comparison to the original BELLE work [3| is the common description of 
both the 7i7i and KK channels and the inclusion of a background amplitude in 
both channels with a phase according to a broad resonance. A more complete 
analysis should provide a fit of the two-dimensional Dalitz plot density taking 
into account the interferences of resonances in the crossed channels. Here we 
are only interested in the low mass region below 1700 MeV. In this region, 
we may neglect the small contribution from crossed channel resonances in 
K + 7t + tt~ (K*); in K + K~ there seems to be some additional contribution 

4 We thank A. Garmash for clarification of this point. 
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below the low mass enhancement which could come from charmed meson 
sources in the crossed channel or some other background. 

In any case, we consider only single channel amplitudes (with interactions 
in 7r + 7r~ or K + K~) which we describe as superposition of the three resonance 
contributions. The spectrum in the pair mass m is then given by: 

— = |ci| 2 gp|T 9b + c 2 T h S gb + c 3 T f ,S gb \ 2 (13) 

TTi r 

T a = — 2 • Vn i n i \\ > a = 9b,fo,fo (14) 

m 2 a -m 2 - im a T a {\ + G a (m)) 

S gb = e 2iS »o; T gb =\T gb \e iS »*. (15) 

The superposition of a narrow resonance with a broad background we 
describe by the rotation of the resonance term in the complex plane with the 
background S matrix, S = e 2lSbg , which is consistent with unitarity in the 
elastic region as is well known; here we also apply it in the inelastic region. In 
(|13|) q is the hadron momentum in the h + h~ restframe and p the momentum 
of the h + h~ pair in the B restframe. 

The resonances included here cannot be given by a simple Breit Wigner 
form in the considered mass range because of the distortion by inelastic KK, 
rjr) and 4tc (pp) thresholds. In our present exploratory study we do not 
attempt a full treatment of unitarity effects of the coupled channel system, 
rather we allow for an energy dependent width which takes into account the 
KK threshold. As we are unable to fix the shape function G a (m) in our 
fitting of projected densities we choose, for definitness, a form which ties the 
amplitude at energies > 1.6 GeV back to the original resonance circle and is 
otherwise adjusted to give a reasonable representation of the mass spectra. 
For / (980) and the glueball (gb) we chose 

Kirn : G a (m) = Q{m — 2niK)s- K x exp(— (m/mi) 11 ) (16) 
KKK : G a {m) = 6(m - 2m K )e K x exp(-(m/m 2 ) 1L5 ) (17) 

with x = yjl — Am 2 K /m? and parameters e n = 0.9, Ek = 2.5, mi=1.28 GeV, 
m2=1.37 GeV. For / (1500), further away from the important thresholds, we 
put G a (m) = 0. The resonance parameters and production amplitudes are 
taken as in Table El Concerning the glueball parameters we note that the 
elastic tttt scattering amplitude in the intermediate energy region around 1 
GeV can be represented by a superposition of T-matrix poles from / (980) 
and a broad state with mass near 1000 MeV and a width of 500-1000 MeV 
as is known since long |33|34j . More recent analyses prefer a slightly higher 
mass in the region 1200 - 1600 MeV for this broad state by fitting to a larger 
number of channels |42j . 
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Table 5: Parameters for resonances and their production amplitudes used in 
the description of the tttx and KK mass spectra in Eq. ()13j) . 



state 


mass [GeV] 


width [GeV] T(n+ir-) 


T(K + K-) 


/o(980) 


0.99 


0.05 c 2 = 5.00 e 0Ain 


c 2 = 2.75 e 0Ain 


/o(1500) 


1.48 


0.10 c 3 = 1.35 e L50 ^ 


c 3 = 0.34 e 75i7r 


gb 


1.15 


0.70 1 


1 



The amplitudes were obtained by adjusting the parameters "by eye" to in- 
vestigate the result from superimposing the 3 resonances with relative phases 
where we kept the branching ratio of / (1500) at the value suggested by the 
PDG (i.e. K + K-/n + 7i- = 3/4 x 0.241) and to describe the data at least 
qualitatively. We do not include the production of the p meson which is 
clearly present [3] . The results with these parameters are shown in Figs. 141 ">l 
in comparison with the BELLE data. 

We note that the relative phase between / (1500) in tttt and KK is 0.757T, 
close to the expected value ti for a near octet resonance, see Table 03 This 
phase difference causes the peak and dip in K + K~ and 7t + tt~ respectively 
which we consider as an important further confirmation of the near octet 
nature of this meson which does not allow a sizable glueball admixture (the 
glueball would contribute with the same sign to both particle pairs, see also 

In our parametrisation the phases for /o(980) — > 7r + 7r and /o(980) — > 
K + K~ decays are about the same, but these phases are less well determined 
by the data. Note that the relative phase between both / (980) amplitudes 
is about zero relative to the background in Table |3J The glueball amplitude 
representing the broad background has phase 6 g b ~ tt/2 near the pole position 
around 1000 MeV. Therefore the interference with /o(980) is destructive 
because of the phase factor S ~ e m in ()13)) . This is well known for elastic 
mi scattering where /o(980) produces a dip instead of a peak in the mass 
spectrum. In our present fit the small background doesn't produce a dip but 
it reduces the size of the / (980) peak in Fig. HJ however, the phase is not 
so well determined and a constructive interference cannot be excluded. This 
requires a study of the lineshape with higher accuracy. Such a phenomenon 
is met also in the decay J/ip — > (finn [72J where the interference of /o(980) 
with the background is destructive. This is derived from the asymmetric 
shape of / (980). 
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Our description also reproduces approximately the small enhancement 
near 1300 MeV which is generated by the interference of all three states. 
On the other hand, it has been interpreted by BELLE as due to /o(1370). 
Whether /o(1370) really exists as a particle is still controversial. The PDG 
does not confirm any branching ratio. If the measurement of the decay ratio 
B(/o(1370) -> K K)/B(f (1370) -> tttt) = 0.46 ± 0.15 ± 0.11 [71J is taken 
for granted, one should see in B decays in both channels an effect of size 
comparable to /o(1500). The K + K~ spectrum does not provide immediate 
evidence for such a situation. 

Finally, we determine the decay rates for the resonances studied and 




Figure 4: 7r + 7r~ mass spectrum in B-decays as measured by BELLE jHj in 
comparison with a model amplitude \T\ 2 of the coherent superposition of 
/o(980), /o(1500) and a glueball (gb) forming the broad background. Also 
shown are the individual resonance terms |Tr| 2 . The background (gb) in this 
fit interferes destructively with both / (980) and / (1500) consistent with 
the model (Eqs. (gj) and (H). 
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M ( K + K 



Figure 5: K + K~ mass spectrum in 5-decays as measured by BELLE |3] in 
comparison with the model amplitude, see Fig. (@J). Here /o(1500) interferes 
constructively with the background. 

compare with our expectations from Table EJ From our analysis we obtain 
the number of events for each state (the area under the curves for our three 
states |T a | 2 in Figs. E]andEJ). Relating the total event numbers in the two 
channels with the corresponding branching ratios for charmless decays we 
arrive at the numbers in Table El In order to obtain the total rates we 
added for / (980) the observed decays in im and KK. In case of / (1500) a 
complete sequence of branching ratios is not available from the PDG. If we 
use the ratios determined by the PDG j23] and estimate the ratio B(fo —* 
47r)/£>(/o — * 27r) pa 1.7 from their listing we obtain the following estimate for 
the branching ratios 

/o(1500) : 4vr (54.2%), 2vr (32.0%), KK (7.8%), rjrj (4.2%), ip{ (1.8%). 

From the 2ir decay we obtain the branching ratios of Table |H1 
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Table 6: Decay branching ratios for resonant contributions in B + — ► 
K + 7T + 7t~ and B + — > K + K + K~ for masses m < 1.72 GeV from our analysis 
using the branching ratio equivalents of 0.05417 x 10~ 6 and 0.03657 x 10 -6 
per event (as concluded from for both final states respectively. 



channel gb [10" 6 ] / (980) [lO^ 6 ] / (1500) [10" 6 ] 



7T + 7T 


7.4 


19.3 


4.3 


7T7T 


11.1 


29.0 


6.5 


K + K- 


10.1 


2.1 


0.70 


KK 


20.2 


4.2 


1.4 


all 




33.2 


20.3 



We note that our result for / (980) — > it + it~ is about twice as large as 
in Table 0] whereas our result for / (1500) — > K + K~ is smaller by the same 
amount compared to Fit Bq / II and much smaller than Fit Bq/I in Table EJ 
Note that in case of /o(1500) we fit the tttc spectrum as well with the KK/rni 
ratio kept fixed (at PDG value) and the large rate of Fit Bq/I seems to be 
excluded. 

If we had taken the fit results by BELLE on / (980) and / (1500) (Fit 
Bq/I I) instead of ours in in Table El we would have obtained the correspond- 
ing rates 17 x 10" 6 and 33 x 10" 6 (instead of 33.2 x 10" 6 and 20.3 x 10" 6 ). 
The different results may come partly from a different treatment of phases, 
especially our moving background phase; we also remark that our analysis 
is using an approximation by working only with projections. We estimate 
that the present numbers for the rates have a model uncertainty by about a 
factor of two. 

Finally, we compare our results with the expectations from Table 01 based 
on symmetry relations for the nonet. We derive the predictions using the 
mixing angle sin cp = l/v3 as in our model [32] and i n Ell- We are using the 
simplifications motivated by the experience with pseudoscalars as outlined in 
the last subsection. We assume real jps and restrict to [3ps = ±1. Then we 
can use the rates for f = / (980) and f' Q = /o(1500) to determine the two 
parameters |pps| 2 and jps- The result for Ippsl 2 may then be compared with 
the measurement of the i^o(1430)7r rate. We allow for both signs of f3ps and 
obtain for each case two solutions from the predictions for the rates following 
from Table El (amplitudes in brackets) and these results are presented in 
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Table 



Table 7: Parameters |pps| 2 (in units of 10~ 6 ) and 7p^ in the decay of B 
into scalar particles using the symmetry relations in Table H3 and the total 
branching ratios for a , f = / (980) and f' Q = /o(1500) from Table E3 There 
are two solutions for each value of ftps] further predictions of branching 
ratios (in units of 10~ 6 ) using for if* (890): •yvs — Jps an d Pvs — ~Pps an d 
Pvs/Pps — 0.437; Sol. 2 excluded, see text. 



Sol. 


I3ps 


\Pps\ 2 


IPS 


K°a+ 


K* + f 


K* + f 


k*°4 


1 


+1 


32.0 


-1.38 


32.0 


13.2 


1.8 


14.0 


2 


+1 


234. 


-0.51 










3 


-1 


11.0 


1.33 


11.0 


55.4 


2.8 


4.8 


4 


-1 


42.1 


-0.30 


42.1 


9.9 


0.6 


18.4 



With \/3ps\ = 1 we can compare these results for |pps| 2 directly to the 
rate for B + — > ifQ°(1430)7r + according to Table El Correcting the numbers 
in Table 0] for neutral decays one finds for Kq°7t + the rates about 38 x 10 -6 
or 9xl0~ 6 . Then we can exclude our Sol. 2 in Table whereas Sol. 1 and 
4 are near Fit Co /I in Table |U and Sol. 3 near Fit Co/ II. A more detailed 
investigation of the Dalitz plot densities could possibly distinguish between 
Solutions C / 1 and Co/ II- Another test is possible through measurement of 
the K°a~Q decay which is given directly by |pp,s| 2 . 

We note that Sol. 4 in Table [3 has similar parameters to the pseudoscalar 
amplitudes (3pp = 1 and 7pp = 0.44, except for their opposite sign. On the 
other hand, Sol. 1 and 3 have a three times larger gluonic amplitude which 
would not be theoretically expected. If we had taken the BELLE Fit results 
instead of ours in Table El we had obtained the corresponding numbers to 
Sol. 4 as \p PS \ 2 = 30.5 x 10~ 6 and 7 P5 = -0.20. 

In any case, the results obtained are not in contradiction with the sug- 
gested scalar nonet of heavier particles along Route A. Actually, our specific 
approach [^2] with nonet ©,© and broad glueball explains naturally the 
new observations from B decays. Further results concern the predictions 
for production of do and of scalars together with vector mesons, especially 
fo, /q. These predictions are shown also in Tabled for the 3 remaining so- 
lutions. Again we make the simplified choices motivated by the experience 
with the pseudoscalars, as outlined in the previous subsection, and assume 
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Ivs — Ips an d Pvs — —ftps (the opposite sign for /3 VS comes from the 
spin factor (— 1) L ), also we take Pvs/Pps — 0.661. Then we expect a much 
smaller rate for K* /o(1500) than for i^*/o(980). The rate K*°a^ determines 
directly the normalization |py,s| 2 - These measurements will provide further 
important tests of our picture. 

3.3 Comparison with light nonet scheme 

In this approach / (980) falls into the same nonet with ao(980) and also with 
a (600) and k(850) whose existence are still under debate. If this nonet is 
built from qq states with mixing as in (jllj) we can apply the same discussion 
as before and explore the symmetry relations of Table El In case of a qqqq 
model additional degrees of freedom may come in which we do not consider 
here. As emphazised in the beginning of this section, there are not yet any 
definitive observations of these light scalars in B decays, as they exist for 
the heavier scalars, nor are there any limits for branching ratios. Therefore 
we only indicate some possible scenarios to stress the potential of i?-decay 
studies also for these light scalars. 

In this case only the rate for B + — > K + /o(980) is available which we take 
from ()1()J) . A natural choice is the mixing angle <ps — which, according to 
(fTTj) . corresponds to / (980) being a pure strange and a a pure non-strange 
state. The decay amplitudes follow from Table E3 and are listed for a few 
decay channels in Table |H1 

The simplest choice is to assume a production without gluonic processes, 
i.e. 7P5 = (Scenario I) We assume again \flps\ = 1, as in the previous 
discussions of pseudoscalar and scalar sectors. Then, with the fo rate as 
input we can predict the other states of the multiplet. One finds, that in this 
case the Ka rate is 1/2 of the Kf rate and equal to the Kk rate. This looks 
rather large for a in view of first results [Sj. 

We therefore also consider another Scenario II where we introduce a glu- 
onic coupling of fo such as to cancel the Ka decay. In this case the sign 
of ftps matters and we obtain for Kk either 60 x 10~ 6 or 6.7 x 1CT 6 where 
the first choice can presumably be excluded from the data [3]. Furthermore, 
we present some predictions for decays with vector meson K*. Again, the 
measurement of the ao rate would be very useful as it fixes the normaliza- 
tions for the PS and VS multiplets. We conclude that a measurement of 
the rates for a and k will be important for the discussion of their existence, 
their classification and their production mechanism. 
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Table 8: Rates for some channels with scalars assuming / (980) in a nonet 
with <t(600) and k(850) (mixing angle ips = 0) and amplitudes (in units of 
the penguin amplitudes pps, Pvs) for two scenarios with different gluonic 
component 7ps, using Pps = ±1. For fP-decays we take Pvs = ~Pps, 
Ivs = Ips and pps/pvs = 0.661. The K + k° and K°clq rates equal |pps| 2 . 



B+ -> 


K + f 




K°4 




K* + f 


K* + a 


amplitude 


Pps + Ips 


^(1 + 2 7 ps) 


1 


Pps 


Pvs + Ivs 




Scenario I : 


Ips = 












rate [10~ 6 ] 


15. 


7.5 


15. 


15. 


6.5 


3.2 


Scenario II : 


Ips = -0.5 


Pps = +1 










rate [1(T 6 ] 


15. 


0. 


60. 


60. 


59.0 


0. 






Pps = -1 










rate [10~ 6 ] 


15. 


0. 


6.7 


6.7 


0.7 
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3.4 Total rate for gluonic decays 

Next we compare the gluonic production rates for f , f' Q and rj, rf with the 
total rate b — > sg in (JIJ). CLEO jTTj has measured the inclusive non-charm 
decay B[B — > ry'+X) = (6.2^|g) x 10~ 4 , where the signal refers to the region 
2.0 < p v > < 2.7 GeV of the ?/ momentum. Identifying the non-charm rate 
with X s according to the SM and adding the exclusive rj'K rate we obtain 
the inclusive rate B(B — > r( + X s ) ~ 7.0 x 10 -4 , so the ratio of the total 
inclusive rj'X s over the exclusive rj'K rate is ^/(incl/excl) ~ 9. We take the 
gluonic part S^psPpsI 2 as in (JHJ whereas for scalars we find from Tabled 
the corresponding result in the range (11 . . . 182) x 10~ 6 with preference for 
the lowest value. Then we find for the fully inclusive contribution of these 
decays after multiplication with iiy(incl/excl) 

B(B - V , V ', f Q , f' Q ) | g i uonic ~ (0.2 . . . 2) x 1(T 4 (18) 

with preference for the lower value. Hence, these decays cannot contribute 
more than a small fraction of the expected b — > sg rate of 5 x 10 -3 [2*U] . 
We will argue next that glueball production does provide the dominant part 
of the b —> sg decay originating from a genuine hard gluon composing the 
associated local dimension 5 operator. 
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4 Glueball production in B decays 



Besides the observation of the strong /o(980) signal there is another inter- 
esting clear feature in the B decays: the presence of a broad low mass tttt 
and KK enhancement in B + — > K + -k + -k^ and B + — > K + K"K + with spin 
J = 0, observed by the BELLE collaboration |1I2I3| . 

An important signal for the sizable S wave background is the interfer- 
ence with known resonances. The appearence of the /o(1500) peak in K + K~ 
on one hand and its disappearence in 7t + tt~ on the other hand can be nat- 
urally explained by its constructive and destructive interference with this 
background as discussed in the previous section. 

This interference is quite pronounced, as in elastic nir scattering, where in 
contrast /o(980) and /o(1500) appear both as dips in the broad background 
which we interpreted in our earlier study as destructive interference with 
the broad glueball ("red dragon"). In the low energy region below 1 GeV 
we assume that the tctt amplitudes are moving as in elastic tttc scattering, 
at higher masses the movement of the phase is more difficult to predict be- 
cause of the inelastic channels. The study of interferences in the Dalitz plot, 
especially the crossing regions of the known resonances - charmed and un- 
charmed - with the background and among themselves could eventually show 
whether the movement of the background phases in tttt and KK channels is 
consistent with an inelastic broad resonance. An open question here is also 
the relevance of /o(1370) in these channels. We argued in Ref. [22] that the 
Breit-Wigner phase motion of an assumed resonance with moderate width 
at this energy has not been demonstrated clearly enough to require a defi- 
nite resonance. In a recent analysis of the /o(1370) — > An channel [73] (as 
reported in 74,) it was found that in the considered mass range the phase 
motion was lower than expected from the Breit Wigner formula so that these 
data do not support the interpretation of /o(1370) as "normal resonance". 

Besides the prediction of a moving phase the glueball hypothesis has also 
consequences for the decay fractions into different particles. In order to relate 
different channels and to obtain an estimate of the total glueball production 
rate we consider the following decay scheme. The glueball decays first into 
the various qq pairs with equal amplitude (possibly also into a pair of gluons) 

gb — > uu + dd + ss {+gg)- (19) 

The gg pair could hadronize into two secondary glueballs after creation of 
another gluon pair. Our ++ glueball at energies above 1 GeV could decay 
as gb — ► aa where <r(600) is considered as the low mass part of the same ++ 
glueball. Then the main decay process would be gb — > An . In the following 
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we do not consider these decays further here. Each of the qq pairs in (fT§j) 
recombines with a newly created pair uu, dd or ss where ss is produced with 
amplitude S (\S\ < 1). In this way the 2-body channels gb — > qq' + qq' are 
opened, at low energies just pairs of pseudoscalars. They are produced with 
probabilities 

7T+7T- ttV K + K- K°K° rm rj'rj' 

2 12 2 1 1 

(20) 

, n 2 1 I|l + 5| 2 ±|2 + S| 2 ±|1 + 2S| 2 

The first row corresponds to U(3f symmetry (S = 1), the second row to 

arbitrary S (for numerical estimates we take S = 0.8); rj, rj' mixing is assumed 

as above. With increasing glueball mass the qq pairs can decay also into pairs 

of vector mesons or of other states. The total rates in (J2Uj) are assumed to 

remain unaltered but in general u tttc" is meant to include pp as well above 

the respective threshold of about 1300 MeV. 

We consider first the mass region 1.0-1.7 GeV. In this region the pseu- 
doscalars alone saturate the U KK" rate in (|2()jl as K*K is forbidden by parity 
and K*K is kinematically suppressed. Another possible decay is 7/77, contri- 
butions from higher mass isoscalars (uu) are only possible at the upper edge 
of the considered mass interval. The decay rfrf is kinematically forbidden. 

In the region 1.1 - 1.3 GeV, away from the narrow resonances and from 
the major inelastic thresholds, we may compare the rates into the 7r + 7r~ 
and K + K~ channels with the theoretical expectations ([20)1 . From Figs. 
0] and El we obtain dR/dM(ir + 'n-~) ~ 20 events/50 MeV in the mean and 
dR/dM(K + K-) ~ 30 events/50 MeV corresponding to dR/dM ~ lx 
10 _6 /50 MeV in both cases using the £>-equivalents in Table El The glueball 
fractions in our fits appear at the level of about 1/2 of the experimental 
data in the mean in both cases in this mass range. Therefore the glueball 
decay rates in 7r + 7r~ and K + K~ are approximately equal, within about 20%, 
consistent with the expectations 2 : 1.6 (for S = 0.8) from (j20j) . 

Above 1600 MeV the two spectra look quite different: whereas the KK 
spectrum decreases slowly, the 7T7T spectrum stays at a low level. This may 
be explained by assuming the dominant decay of "7T7t" in (|2Uj) proceeds into 
4tc (e.g. pp) states in this mass range. 

Next we estimate the total glueball rate. Here we start from our fit 
result for gb — > KK in Table El refering to the mass interval 1.0... 1.72 
GeV. Including the contribution from higher masses of about 55% we obtain 
B(B + — > K + gb; gb — > KK) m 31 x 10~ 6 . This yields an estimate of the lower 
limit for exclusive and inclusive branching ratios (correcting for other decay 
modes without gb gb by a factor 2.2 from (j2U|) using S = 0.8 by neglecting 
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7T7T decays below 1 GeV and r/r/') 



B(B 



+ 



gb(0 ++ ) + K + ) > 70 x 10 



-6 



(21) 



B{B 



+ 



gb(0 ++ )+X s ) > 0.6x10 



,-3 



(22) 



here we used again the factor 9 from Krj' to estimate the fully inclusive rate. 

By adding the gluonic pseudoscalar and scalar meson contributions from 
()18|) we estimate the lower limit (neglecting again the gg decay mode in (|T9"j) 
and choosing Sol. 4 in Table |7J) for the total production of observed gluonic 
scalar and pseudoscalar mesons as 



B{B + -> gb(0 ++ ) +f + ti + r ] + r] , + X s ) | gluonic > 0.8 x 10" 3 . (23) 



This lower limit for gluonic production amounts to about 1/2 of the leading 
order result for the process b — > sg in ((TJ) and about 1/6 of the full rate ob- 
tained in NLO, so it represents already a sizable fraction of the theoretically 
derived value. This result supports the expectation that the b — > sg rate will 
be saturated if a few similar processes with other J PC quantum numbers, 
especially the 0~ + and 2 ++ glueballs are included. 

5 Conclusions 

The main motivation for our investigation is the search for the lightest glue- 
ball which is expected with J PC = ++ quantum numbers. To this end it 
is important to obtain a full understanding of the light scalar sector, i.e. to 
establish the lightest scalar nonet and its intrinsic mixing. An important 
role is played by /o(980) which in different classification schemes is either 
the lightest or the heaviest particle in this nonet (other options are being 
discussed as well). 

In this paper we investigate the potential of 5-decays which have the 
partial decay mode b — > sg in the search of gluonic objects. This interest 
has been triggered by the undisputable observation of B — > Kf o(980) with a 
large rate comparable to the one for pseudoscalar particles. This could imply 
a large gluonic affinity of this meson or a flavour singlet nature similar to 
the r]'. 

For the further investigation of -B-decays we suggest a simple approximate 
scheme for 2-body decay rates based on the dominance of penguin amplitudes 
with an additional gluonic component; for each pair of multiplets there are 3 
amplitudes p, jp, ftp, the present analysis allows for real 7 and ft = ±1. This 
model has been tested first in the sector of decays B — ► PP and B — > VP 
where it corresponds to a simplified version of the previous approach [HI] . 
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Because of the (approximate) flavour symmetry of all qq in the decays b — ► 
sqq all members of the nonet have a common component in the amplitude 
modified by ft, 7 amplitudes. Exploiting this fact by comparing the rates 
with the expectations in Table El should ultimately disclose the identity of 
the members of the lightest nonet associated with /o(980). 

The choice of the nonet with members a (980), / (980), ^(1430), / (1500) 
with / (980) as lightest particle as in [32] and similarly in [37| (with a (1450) 
preferred) can reproduce the observed phenomena concerning also other 
scalar particles not yet as well established as /o(980). In particular, there is 
the remarkable phenomenon of the /o(1500) signal in KK, but apparently 
absent in despite the four times larger branching ratio. In our analysis 
this is explained by the constructive and destructive interference respectively 
with the broad glueball ("background", "red dragon"). Such a behaviour is 
expected from the near octet flavour composition of /q(1500). The negative 
sign between the ss and uu + dd component of /o(1500) has been found before 
in inelastic tttt scattering [22]. If confirmed in the final analysis of the data 
it would seriously restrict the possible glueball admixtures of / (1500) which 
yield contributions of equal sign to tctc and KK decay amplitudes. 

An important role in this classification is played by ao(980) as the lightest 
isovector particle which directly determines the penguin amplitude p and 
therefore the overall normalization within one multiplet. This particle can 
be identified from B —> KKK and B — > KrjTr decays. Knowing this decay 
rate the other parameters in the model for the given decay multiplets can be 
determined more directly. Alternatively, ao(1450) could be the nonet partner 
of f (980). 

The a (980) rate would help in particular in a judgement about the clas- 
sification with / (980) as heaviest particle in the nonet together with a and 
k in 5-decays. At present there is no strong indication for the presence of 
these particles but a more dedicated analysis in the determination of the 
respective decay rates is necessary. 

The study of charmless B decays into the various members of the nonet 
represents a systematic approach to scalar spectroscopy which can be fol- 
lowed in great analogy to the successful phenomenology of the decays into 
pseudoscalars. A particularly interesting test is the comparison of decays PS 
and VS with K and K* which are expected in some cases with quite differ- 
ent rates because of different signs of the parameter ft. Our first attempts 
indicate a different sign of the amplitudes 7, ft in final states with scalar and 
pseudoscalar particles. 

Finally, we come back to the question of gluonic meson production. First 
there is the gluonic component in the production of isoscalar mesons which 
is obtained for 77, 7/ already in previous analyses. We find such components 
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also in the scalar sector for / (980). A more definitive analysis can be done 
given the a rates. 

The presence of a coherent background can be derived from its interfer- 
ence with / (1500). In our glueball interpretation we expect a production 
phase moving slowly with energy according to a Breit Wigner amplitude with 
some modification by inelastic effects. It will be interesting to determine more 
accurately the interference between the background (gb) and / (980) whether 
it is destructive (as in J/tp — > (pmr and in our model) or constructive; this has 
a big influence on the derivation of the important /o(980) decay rate. The 
7T7T / KK branching ratios of the background from our estimate are consistent 
with this glueball hypothesis. Furthermore we predict a sizable rji] and, at 
higher masses, a An decay rate. 

The broad background has also shown up in other gluonic processes like 
double pomeron production and pp annihilation, also with suppressed rate 
in 77 collisions |75|7fij : on the other hand, it has been difficult to see its sign 
in J/ij) — > 7T7T7, KKy (see discussion in [H2])- Recent high statistics results 
by BES [77j on KKy, however, require a broad ++ coherent background 
for a good description of the data as well. We have no definitive explanation 
why B(J/ip — > /o(980)7) is apparently much smaller than B(J/ip — > rfy) 
despite the similar quark structure proposed for /o(980) and rj . A possible 
explanation is a destructive interference with the background, seen in other 
processes as discussed, such that the appearence of / (980) is minimized, also 
a special effect in the exclusive decay could be thought of. Therefore, it would 
be interesting to see whether the similarity of rf and /o(980) is recovered in 
the fragmentation region of gluon jets as discussed in 0. 

It remains an interesting question how the b — > sg decay is realized by 
hadronic final states. The large rate for the ++ glueball we obtain within 
our approach suggests the intriguing possibility that it could be saturated by 
gluonic mesons. In the next step it will be interesting to search for the CT + 
glueball which could decay into rj-n-K and KKtt. 
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